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Abstract 

Historically, there have been many attempts to define the correct algebra for modeling the 
properties of three dimensional physical space, such as Descartes' system of Cartesian coordinates 
in 1637, the quaternions of Hamilton representing rotation in three-space that built on the Argand 
diagram for two-space, and Gibbs' vector calculus employing the dot and cross products. We 
illustrate however, that Clifford's geometric algebra developed in 1873, but largely overlooked 
by the science community, provides the simplest and most natural algebra for three-space and 
hence has general applicability to all fields of science and engineering. To support this thesis, 
we firstly show how geometric algebra naturally produces all the properties of complex numbers 
and quaternions and the vector cross product in a single formalism, whilst still maintaining a 
strictly real field and secondly we show in two specific cases how it simplifies analysis in regards 
to electromagnetism and Dirac's equation of quantum mechanics. This approach thus has the 
immediate advantage of removing complex fields from analysis, so that algebraic entities have a 
geometric expression in real space. As an example, we show how quadratic equations now can be 
given two real solutions using GA. This viewpoint also has something interesting to say about the 
concept of number itself, because numbers in GA can encapsulate scalars, lines, areas and volumes, 
into a single entity, which can be multiplied and divided, just like normal numbers. This simple 
and elegant formalism would also seem very appropriate for students learning vectors, algebra and 
geometry, giving a more natural and intuitive understanding of the properties of three-space. 
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I. INTRODUCTION 



Einstein once stated, 'Everything should be made as simple as possible, but not one bit 
simpler', and in this paper we ask the question: 'What is the simplest and most natural 
algebraic representation of three space?' 

It is readily accepted from experience and scientific experiment that we have three spa- 
tial dimensions or independent directions of movement within which to describe events, and 
hence Cartesian coordinates seem eminently appropriate for the task of providing a frame- 
work for three-space, where we typically label the three orthogonal directions as x, y, z. We 
ignore in the first instance the presence of curved space, or other minute curled up dimen- 
sions. For three-space, though, as well as positional coordinates, we also need to be able 
to represent orientation or rotations at each point in the space. In two-space the correct 
algebra is defined by complex numbers and for three space, the algebra for rotations is given 
by Hamilton's quaternions. Quaternions are seen as somewhat unusual due to their non- 
commuting properties, however this reflects the observed non-commutivity of rotations in 
three-space. This can very easily be demonstrated by rotating any object, with a 90° and 
180° orthogonal rotation, where one will then find the relation that RisoRgo = — -Rgo-Riso- 
Hence, in order to form a unified algebra of three-space we need to integrate the complex 
numbers and quaternions within the framework of Cartesian coordinates. This was achieved 
by Clifford in 1873, who named his system. Geometric Algebra, or GA. 

As well as deriving Clifford's system, we also illustrate, how the formalism indeed re- 
places, complex numbers, quaternions, the vector cross product, rotation matrices, the Pauli 
matrices and the anti-symmetric field tensor as well as several other formalisms. Clearly 
this simplification should have general applicability to many fields of science, physics, and 
engineering as well as being advantageous in an educational setting. 

A. Historical development 

Around 547 BCE, Thales expressed his belief that every event on the earth had natural 
rather than mythological causes. For example he suggested that earthquakes were not 
caused by the machinations of the gods on Mount Olympus, but proposed instead that the 
earth floated on water and was being buffeted by large waves [Ij. Considering the limited 
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knowledge of the time, this showed very good insight, as we would perhaps refine this today 
by saying that the earth fioats on molten magma and is buffeted by tectonic plates. For 
this reason Thales is sometimes identified as the first scientist. Pythagoras extended this 
notion of logical causation behind events, postulating that numbers and their relationships, 
or mathematics, underly all things. This idea was then masterfully applied by Euclid to 
geometry, deriving a multitude of geometrical theorems from a few simple axioms. Euclid's 
books became standard texts for the next 2000 years, only disappearing from school curricula 
in the 20th century! 

The next major development came in 1637 when Descartes unified algebra and geom- 
etry onto the Cartesian plane, so that the three independent spatial freedoms, could be 
unambiguously represented by a single set (vector) holding the three coordinates. This 
achievement is called by John Stuart Mill, 'the greatest single step ever made in the exact 
sciences' [2]. It has been extensively debated by historians, however, why there was such 
a slow down in the progress of science and mathematics for over 1500 years following the 
Greek explosion. Various suggestions have been provided to answer this, such as the Ro- 
man empire suppressing dissent and not sponsoring the arts, the rise of the roman catholic 
church, the ready availability of slaves obviating the need for work efficiencies and the need 
to educate them and so on [3J. However Hestenes has put forward another possibility, that 
the algebraic and numerical system used by the Greeks, had inherent limitations, which were 
roadblocks to further progress [1]. 

For example the length along the diagonal of a unit square, we know today as y^, being 
an irrational number, did not exist in the Greek numeric system, which was based solely on 
integers and their ratios. Another hindrance would have been the roman numerals which 
made numeric manipulation difficult. Also the algebraic expression x^, or x cubed, referred 
to a physical cube and hence x^ was unacceptable on idealogical grounds. However with the 
arrival of Hindu- Arabic numbers in about 1000 AD into Europe, which included a zero, thus 
allowing positional representation for numbers in base 10, together with the acceptance of 
negative numbers in 1545 AD, allowing for the concept of a complete number line, paved the 
way for Descartes to revolutionize the Greek system, by proposing a union of algebra and 
geometry using Cartesian coordinates. He stated 'Just as arithmetic consists of only four or 
five operations, namely, addition, subtraction, multiplication, division and the extraction of 
roots, which may be considered a kind of division, so in geometry, to find required lines it is 
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merely necessary to add or subtract lines.' Descartes thus postulated an equivalence between 
line segments and numbers, something the Greeks were not prepared to do. Effectively he 
generalized the concept of number to include directed line segments. We will show in fact, 
that the concept of number can be further extended beyond line segments, to include planes 
and volumes. 

This simple approach of representing two-space and three-space by Descartes, however, 
becomes confused with the development of the Argand diagram, which also represents the 
plane, but consisting of one real and one imaginary coordinate, and to make matters worse, 
Hamilton in 1843 generalized complex numbers to three space, naming them quaternions. 
This confused state of affairs, on exactly how to represent three-space coordinates and rota- 
tions, was finally resolved by William Clifford in 1873, who combined Cartesian coordinates 
and the rotational algebra of both complex numbers and quaternions into single unified 
algebraic framework. 

However Clifford also achieved a fulfillment of Descartes' original vision of a vector being 
able to be manipulated in the same way as normal numbers, by deriving the appropri- 
ate multiplication and division operations for vectors, which meant that geometric objects 
such as vectors could now treated with normal algebraic operations. Cliffords mathemati- 
cal system 'should have gone on to dominate mathematical physics' [5J, but, Clifford died 
young, at the age of just 33 and vector calculus was heavily promoted by Gibbs and rapidly 
became popular, eclipsing Clifford's work, which in comparison appeared strange with its 
non-commuting variables. 

Hestenes has championed GA for over forty years and the development of GA is now 
expanding rapidly, with research in black holes 0, quantum field theory [7], quantum tun- 
neling |8], quantum computing [9], general relativity and cosmology [10], beam dynamics 
and buckling [11], computer vision fT2] and EPR-Bell experiments |13j . 

B. Clifford's definition of three-space 

How did Clifford solve the problem of the correct representation for three-space, inte- 
grating the rotational algebra of complex numbers and quaternions? In order to describe 
Clifford's approach, we firstly define the three Cartesian coordinate directions by ei,e2 and 
63 as shown in Fig [T| as proposed by Descartes. We specify the orthonormality of these 
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three vectors with the expression Cj ■ Cj = Sij, where the conventional vector dot product is 
adopted without alteration in GA, however Clifford then defined each algebraic element Cj 
as anticommuting, such that CiCj = —ejCi. The bivectors 6263, 6163 and 6263 square to minus 
one, for example, (6263)^ = 62636263 = —62626363 = —1, where = 1, and the bivectors now 
become isomorphic to the quaternions of Hamilton. The trivector l = 616263 also squares to 
minus one and commutes with the three orthonormal basis vectors and hence is isomorphic 
to the imaginary number Note that the trivector l has identical properties, but is not 

the same object as the imaginary number a/— 1, normally represented with the symbol i or 
j, so this is indicated by giving it a slightly different symbol l, the Greek iota symbol. 

We can also find CiCj = teijkek, which we call the dual representation, which we see 
creates the Pauli algebra, and hence we can use Clifford's basis vectors to replace the Pauli 
matrices, as they also square to one and anticommute. 




FIG. 1: Clifford's model for three space. 

Thus far, using Clifford's GA, we have already superseded, complex numbers, quaternions 
and the Pauli algebra, into a completely real three-space algebra. We can now also explain, 
the somewhat mysterious imaginary number i = a/— 1, as simply an algebraic stand in for 
two-space rotations, or the trivector l, also an algebraic object that squares to minus one. 
Also, quaternions, ironically, usually seen by students as some frightening beast locked in 
the mathematics basement, should be seen as the most natural of operations, being simply 
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three-space rotations represented by the bivectors iCi, 162, te^ in GA. 
C. Derivation of the Clifford product 

Chfford also defined the geometric product, using simply the distributive law for expand- 
ing brackets. Hence the product of two vectors becomes 

uv (1) 

= (eitii + 62^2 + 63^3) (eiVl + 62^2 + e3f3) 

= UiVi + U2V2 + M3^^3 + iu2V3 - ^^2^3)6263 + (UiV^ - M3f 1)6163 + {U1V2 - ViU2)eie2 

= UiVi + U2V2 + M3^^3 + {{U2V1, - V2Us)ei + {UiVs - ^3^^1)62 + iuiV2 - ViU2)e3) 

= U ■ V + iU X V, 

= u ■ V + u A V, 



using 6^ = 1 and CiCj = —CjCi, which thus provides a definition of the dot and outer 
products in component form. Hence Clifford's vector product naturally integrates with 
common algebraic operations, being simply the algebraic product of two brackets. We can 
also see that for the case of a vector times itself, that the cross product will be zero and 
hence the square of a vector = v ■ v, becomes a scalar quantity. 

Note that the expression u-v + 6u x v, is in the form of a complex- like number, and in the 
same way as we prefer to leave a normal complex number a + ib a.s a single unit rather than 
separating it into real and imaginary parts, it is also preferable to leave u ■ v + iu x v as a 
whole rather than splitting it into the separate dot and cross products. This, in fact, appears 
to be a foundational error of vector analysis, as the separate dot and cross products, now fail 
to integrate with normal algebraic operations, making it more difficult to carry out algebraic 
manipulation, and hence, also harder for students to master. However if we leave the Clifford 
vector product intact one finds that GA can naturally be employed to represent quantum 
mechanics in real three-space coordinates without the use of Pauli or Dirac matrices |14j . 

The cross product as defined by Gibbs also has a problem, in that it only applies in three 
dimensional space, because in four dimensions there is an infinity of perpendicular vectors, 
and so it is preferable to adopt the wedge product in its stead, which applies generally in 
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any number of dimensions, as does the dot product. Historically, as Gibbs' system of vectors 
became more popular in physics and in various other fields, new scientific discoveries such 
as quantum mechanics and relativity meant that vector analysis needed to be supplemented 
by many other mathematical techniques such as: tensors, spinors, matrix algebra, Hilbert 
spaces, differential forms etc. and as noted in [IS], 'The result is a bewildering plethora of 
mathematical techniques which require much learning and teaching, which tend to fragment 
the subject and which embody wasteful overlaps and requirements of translation.' 



D. Maxwell's equations in GA 

Electromagnetism is one of the foundational theories of physics and we find significant 
simplifications if GA is adopted. Maxwell's equations were first published in 1865 [TB] math- 
ematically describing then recently discovered electromagnetic phenomena. His equations 
were written for three-space, requiring 12 equations in 12 unknowns. These equations were 
later rewritten by Heaviside and Gibbs, in the newly developed formalism of dot and cross 
products, which reduced them to the four equations now seen in most modern textbooks 
[T7] and shown below in S.I. units 

V-E = -, (Gauss' law); (2) 



V X E + dtB = 0, (Faraday's law); 

X B ^i9tE = fiQj, (Ampere's law); 

V ■ B = 0, (Gauss' law of magnetism). 



where E,B,J are conventional vector fields, with E the electric field strength and B the 
magnetic field strength and V is the three-vector gradient given by 

„ d d d . , 

V = ei— + e2T^ + e37r. (3 
ox oy oz 



and 61,62,63 are the three Euclidean space orthonormal vectors, with 6i • Cj = Sij, where 
5ij is the Kronecker-delta symbol. Maxwell's four equations given in Eq. ^ along with the 
Lorentz force law 

K = g(E + vxB), (4) 
8 



completely summarize classical electrodynamics [T7]. It can be seen that we have adopted 
the convention that a three-space vector v = viCi + ^262 + ^363, is identified in bold, in 
agreement with ISO 80000-2 international standard for mathematical notation. 

However inspecting the form of the geometric product for uv = u ■ v + 6U x v, we can 
see that these equations can be fairly easily combined. If we firstly multiply the second and 
fourth equations by l and insert l into Eq. ([s]) as shown 

V-E = ^, (Gauss' law); (5) 

iV X E + dttB = 0, (Faraday's law); 

iV X + -^dtE = — /ioJ, (Ampere's law); 

V ■ iB = 0, (Gauss' law of magnetism). 



We now see that using Eq. Q that we can add the first and second and the third and fourth 
equations to find 

VE + dtiB = ^; (6) 



where VE will now represent the geometric product of two vectors, that is VE = V ■ E + 
iV X E. However, these two remaining equations can now be added to produce 

(^9t + V)(E + icB) = ^-c/ioJ. (7) 

If we define F = E -|- tcB and the four-gradient □ = + V, with the source J = ^ — c/iqJ, 
we find 

OF = J. (8) 

Thus unifying the dot and cross products, naturally reduces Maxwell's four equations down 
into a single equation, demonstrating the power of Clifford's approach. We can see that the 
B field, is written as a pseudovector iB, as part of the field F. This is appropriate, because 
under a refiection, the current in a loop of wire will fiow in the opposite direction, thus 
inverting the B field, whereas an E field is simply refiected about the mirror plane. This 
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different behavior of the E and B fields is obvious from the GA formahsm but hidden in the 
tensor and Gibbs' vector formahsm. 

Also if we wished to describe Maxwell's original four equations in English we would find 
it fairly difficult. For example, we could describe Maxwell's equations as follows: equate the 
divergence of the E field to the charge present, simultaneously equating the Curl of E as 
the negative rate of change with time of the B field, simultaneously, and so on for the two 
remaining Maxwell equations shown in Eq. (|2]). However with GA we can simply say that 
the four-gradient of the field F is proportional to the amount of charge and current present. 



E. The versatile multivector 



As previously shown, the geometric product, is simply an algebraic expansion of two 
brackets, hence there is now no reason why we cannot combine vectors with scalars, bivectors 
or trivectors. 

For example, adding all these components together, we produce a multivector, which can 
be written 

M = a + V + tw + 66, (9) 

which shows in sequence, scalar, vector, bivector and trivector terms. This general three- 
space multivector can be used to represent a diverse range of mathematical objects. Firstly, 
we can obviously map into scalars (a), vectors (v) and complex-type numbers (a + ih). The 
quaternions i, j, k, defined by i^ = = k^ = ijk = — 1, map into the multivector iw, that is, 
component wise, i — )■ tei,j — )■ — te2,k — )■ 663. Hence a general quaternion is mapped to GA 
as 

g = a + 6i — cj + (ik-H-a-|- ihei + Lce2 + Lde^ = a + lw. (10) 
This is similar to the Pauli spinor mapping 



a + jas 
-02 + jai 



■ir^ = a + miCi + ^262 + ^-0363 = a + tw, (11) 



also mapping to the even sub algebra of the multivector, which shows the close relationship 
between spinors and quaternions |T8]. 

The electromagnetic field represented by the antisymmetric tensor F'^'^ [T7], maps as 
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O flo + ^1^1 + &2e2 + &3e3 + 60161 + 60262 + ^0363 + 66 = a + v + 6w + 66. (13) 



follows 

F^"' ^F = E + 6cB, (12) 

with the dual tensor G^'^ given in GA by G = lF. We also have other dual spaces in the 
multivector, such as for vectors v o 6W and for spinors a + 6W o 6(6 — 6v). 

The Dirac wave function, used in the Dirac equation, maps to the full multivector as 
follows 

a + jas 
0-2 + jfli 
-63 + jh 
-61 + j62_ 

Other entities can also be mapped to a multivector, such as pseudoscalars (66) and pseudo 
vectors (6w). Rotation matrices are also represented in a multivector as i? = a + 6w, which 
rotates a vector about the w axis by an angle given by tan^ = \w\/a according to 

y' = RyR\ (14) 

where = CqRcq, where the tilde operation reverses the order of products and 60 represents 
the time coordinate, which anticommutes with the space coordinates 6j. Hence the great 
versatility of the three-space multivector is demonstrated, being able to replace a large vari- 
ety of mathematical structures and formalisms as shown and hence unify the mathematical 
language of physics [H]. 



F. Common algebraic operations on a multivector 

Descartes claimed that five common algebraic operations, addition, subtraction, multipli- 
cation, division and the square root, could be applied to his line segments, however this idea 
can be extended to full multivectors, which contain not only line segments (vectors) but also 
areas (bivectors) and volumes (trivectors). Addition and subtraction are simply defined by 
adding hke components, that is, if Mi = ai + vi + 6Wi + 661 and M2 = 02 + V2 + 6W2 + 662, 
then 



Ml + M2 = (ai + 02) + (vi + V2) + 6(wi + W2) + 6(61 + 62), 
and similarly for subtraction. 
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(15) 



1. Multiplication of multivectors 



The multiphcation operation is given by an algebraic product, similar to the algebraic 
product of two vectors, that is 

M1M2 = (ai + vi + twi + ibi){a2 + V2 + tW2 + ^62) (16) 
= (aia2 + vi ■ V2 — wi ■ W2 — &162 + 02V1 + aiV2 — 62W1 — 61 W2 + vi x W2 + wi x V2 
+ 6(a2Wi + aiW2 + 62V1 + 61V2 + f 1 X f 2 - wi X ^2) + ^(0162 + 0-2^ + Vi ■ W2 + Wi ■ V2)) , 

where we have used repeatedly the geometric product defined in Eq. ([T]). 

2. The inverse of a multivector 

For the division operation, we firstly need an inverse operation for a general multivector. 
For the general multivector defined in Eq. (|9|, we find the inverse 



= (a - V - 6w + ih)U + ig)/{f + (17) 

where we have the scalars f = a? — v'^ + — h"^ and (7 = 2(v ■ w — ah). 
The special case of the inverse of a vector is given by 

= -2. (18) 

which can be checked through vv^^ = \^ = ^ = las required. Hence, we can see that 
the inverse of any vector is a vector in the same direction with the reciprocal of the length 
of the original vector. 

3. The square root of a multivector 

For the multivector defined in Eq. ([9]), we find the square root 

M5 = ±(c + X + 6y + id) (/ + Lg) , (19) 

where we find c = rcosh| and d = rsinh|, where x = ^'("2^™) cind y = 2{c^+^) ^^"^ 
f + '^Q = ^ (^V a + Vh + L^/—a + ^/h^ , with h = a? + IP', using the following support 
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formulae 



1 1 



av + 6w aw — 6v 

m = — 7. — TTT , n = 



a- + 62 ' a2 + 62 ■ 



2^-0/ u. -r f 

(21) 



Because multivector multiplication is associative we can now find all the rational powers 
M'^/2*^ where A;, t are integers, which can approximate to any fractional power. We find that 
the square root generally exists, and in some cases there is a double root (ignoring signs). 



G. Interpreting solutions of quadratics using GA 

Imaginary numbers caused a lot of confusion when first discovered as solutions to 
quadratic equations, however with GA the imaginary solutions i = \/— T, are given a 
geometric meaning using L2 = 6162, the bivector of two-space, which also squares to 
minus one. Hence a complex number solution x = a + ib can be written in GA as 
R — re''^^ — r{cos9 + isin^), a rotor in the plane described using just real numbers. In 
two-space we can rotate vectors using v' = Rv. That is if i? = e'-^'^/'^ — — ^62, then if 
V = 62, then v' = 616262 = 61, or a clockwise rotation by ninety degrees. Hence solutions 
of quadratics using complex numbers, imply we are using rotations in the plane, instead of 
simply scaling along the real number line, that is we are increasing the solution space from a 
line to a plane. This also explains why we always have two symmetrical complex solutions, 
if they exist, as they represent ±0 directions for the rotor. Then, if we write the quadratic 
equation as 

-ax^ -bx^c (22) 

and because x now can become a rotor, we can view this as a rotor equation, where the 
action of the LHS, which involves rotating the starting vector is equivalent to a pure scaling 
of a vector by c. that is, if Qi — —ax^ — bx and Q2 — c, then 

Qiv = Q2V, (23) 
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where the real number hne is assumed to be the ei axis, so that v = ei. We can see that 
the rotor operator Q2 simply scales the initial vector by c. Hence we require the action of 
the operator Qi, which allows rotation of the vector in the plane, to return the vector to the 
starting direction and scale it by c. Hence, thinking in terms of adding vectors head to tail, 
then adding the two vectors — ax^v and —bxv, will produce a vector equal to the vector cv. 
That is 

Qiv = (-ai?2 _ bR)v (24) 
= (-arV^^ - r6e'2^)v 
= cv 



where we used the property of exponentials that (e'^^)^ = e^'^^. 

Firstly we see that the e2-component of the vector transformed by Qi must be zero, which 
ensures the vector is parallel to the starting vector, hence we have the equation 

ar^ sin 2^ + r6 sin ^ = 0, (25) 

giving cos 6* — For the ci direction we have 

- ar^ cos 2^ -r6 cos ^ = c (26) 

—ar^ (2 cos^ 6' — 1) — 6r cos = c 
fe' 2 



giving — ^ and hence we have finally, r = and 9 — arccos and we have therefore 
found the two solutions given by 

X = re^'^^, (27) 

which is in fact equivalent to the standard quadratic formula, but only employing real 
numbers. Also we can view the quadratic equation as a vector equation. Multiplying the 
quadratic equation through by ei we find 

— ar^ (cos 6* + 6 sin 6*) ^ci — 6r (cos 6* + t sin 6*) ci = cci (28) 

— ar^wciW — brw — cei 
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where w = cos^ei — sin is an unknown unit vector, with = 1. Multiplying by wei 
we find 

— ar^w — brei = cciwci 

a vector equation, which we can solve by equating ei and 62 components returning the same 
solution as previously. 

In summary the appearance of imaginary numbers in the solution of quadratic equations, 
can be seen as solutions involving rotations in the plane. This idea can therefore be gen- 
eralized to allow rotations in three-space using three-space rotors. This will be equivalent 
to allowing x to be a quaternion, or represented in GA by defining x = e''"^, where u is 
a unit vector. We find that the same solution is obtained as previously of r = and 
9 = arccosTTT^, which shows that the axis of rotation it has no effect, showing we have 
an infinite number of solutions. For the two-space case we were essentially assuming that 
we were rotating about the 63 axis, clearly then, any other axis will also provide a solu- 
tion to the quadratic. We thus have provided a geometric interpretation of the solution of 
quadratics in two-space or three-space, and thus the solutions of the quadratics are always 
real, representing rotations. Other extensions now present themselves, including, letting x 
become a full multivector. 



H. The Dirac equation 

We find using GA, that we can write the free Dirac equation in real three-space as 

UF = —^F*Le3 (29) 

where F = a + E + tB + i6 and F* = cqFcq = a — E + iB — tfe, space conjugation. 

Thus the Dirac equation using three-space GA simply becomes an extension of the free 
Maxwell equation shown in Eq. ([s]), with the addition of a mass (m) term and expansion of 
the field F to a full multivector as shown. This is clearly a dramatic simplification of the 
standard Dirac equation, which is normally considered embedded in four- dimensional space 
time employing 4x4 complex matrices, as shown in Appendix B. Hence we can recognize 
another over-complication in current mathematical formalisms with this example, in that 
the non-commuting property of matrices is used as a substitute for the much simpler non- 
commuting real basis vectors of GA. Also, with the GA form of the Dirac equation, we can 
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see that it describes a multivector field, given by Eq. ([9]), over real three-space, that is, at 
each point in three-space we have a multivector valued field defined. 

II. CONCLUSION 

In this paper we ask the question about the correct algebraic representation of three-space, 
and we conclude that GA provides the simplest and most natural formalism. Geometric 
algebra that uses non-commuting Cartesian coordinates, subsumes the rotational algebra of 
complex numbers and quaternions into a simple algebraic system over the real field. We 
claim that Clifford's GA therefore has general applicability, which we illustrate with the 
two examples of Maxwell's equations and the Dirac equation, both being written as single 
equations in real three dimensional space. This approach thus has the advantage that all 
algebraic entities have an immediate geometric representation in real three-space. 

We argue that the current system of vectors taught to students, which lacks an inverse 
operation with regards to multiplication, should be replaced with the system of vectors 
described in GA, in which multiplication and division of vectors becomes a natural algebraic 
process. We show in Eq. ([T]), how the geometric product described by Clifford is a natural 
extension of the distributive law for brackets incorporating non-commutivity. The non- 
commutivity is required for three-space, simply because rotations in three dimensions are 
non-commuting. 

We also show more generally how Clifford's union of algebra and geometry in three- 
space allows the viewing of a multivector, as shown in Eq. (|9]) as a generalized type of 
number, useful in representing different physical quantities, amenable to the basic operations 
of addition, subtraction, multiplication, division and square root. We showed how a general 
quadratic equation can be solved without recourse to complex numbers, giving the solutions 
geometric meaning as rotations in the plane, which will be beneficial to students learning 
algebra. 

We highlight two foundational errors in the development of the current mathematical 
formalisms. First, the separation of the dot and cross products into two separate products, 
rather than considering them as a single unified product, and second, the use of the Pauli 
and Dirac complex matrices, as a way to implement non-commutivity, in place of the much 
simpler orthonormal real basis vectors of GA. 
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III. APPENDIX 



A. Application to Electromagnetism 

In order to gain a wider comparison between the conventional vector calculus form of 
electromagnetism and the GA form, we list below various expressions in both formalisms 
including the tensor formalism. The GA formalism uses a single variable F to hold the 
electric field strength consisting of a vector for E and a pseudo vector for the B field, 
combining the four Maxwell equations into a single equation. When calculating the energy 
in the fields or the invariants, there is also the advantage of simplicity, in that only a single 
equation is required to produce the energy and momentum or the two invariants, as shown 
in Table |Tj The energy and momentum expression is in a form similar to kinetic energy, that 
is ^eol-Fp. In the other cases listed in the table, we can see in all cases GA provides more 
compact expressions. 

B. Dirac equation 

Schrodinger developed a non-relativistic wave equation in 1925 which was generalized to 
include spin by Pauli. Dirac then sought a relativistic formulation and succeeded in 1928 in 
producing the explicitly covariant Dirac equation 

- ihYdf,^ + mcip = 0, (30) 

which uses the Einstein summation convention, where 
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(31) 

and 

{Y,Y} = 2g^r (32) 

As discussed, if the gamma matrices are considered as basis vectors, we would then have the 
anti commutator {7'^,7''} = 'y'^'-f" + 7*^7^ = 27^^ ■ 7^^, as expected for a set of orthonormal 
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EM theory 


Gibbs' vector calculus 


Tensors 


GA 


Fields 


E,cB 




E + icB 


EM equations 


V - E = f, V X E + 9jB = 

V X B - ^dtB = fioJ, V • B = 




UF = cfioJ 


Charge conservation 


v-J + t = o 


J% = 


□ • J = 


Energy in fields U 


ie(E2 + c2B2) 


e(i707i70 + IF^^F'^'^) 




Poynting vector S/c 


7r7(E X B) 


eF^^F^ 

1 




Invariants 


c2B2 - E2 

^B E 


iF'^^FaB 
\e^/5j5F"^F^^ 




Minkowski Force 


K = 7g(E + V X B),cK° = q-fv-E 


cK>' = -qvuF^'^ 


cK = —qvF 


Potential function A 


E = -VV- B = V X A 


Fal3 = Aa,i3 — Af^^a 


F = nA 


Lorenz gauge 


A = V-A + ^f =0 


A% = Q 


□ • A = 


Potential form 




□M^ = -c^loJ>' 


n'^A = cjjLoJ 



TABLE I: Comparison of mathematical formalisms, showing relative simplicity of mathematical 
expressions in GA. Note: The blanks in several rows in the column for GA are because the relevant 
equation is automatically part of the equation above. For example, there is a blank row below 
the F'^ entry for GA, because the separate scalar and pseudoscalar components of this expression 
produce the two invariants in the vector column. 

basis vectors. Hence, the opinion of many people, that Dirac indirectly rediscovered Clifford's 
geometric algebra with its anti-commuting basis vectors. 

Dirac's complex, four-space equation using 4x4 complex matrices, can be compared with 



the real three-space version, shown in Eq. (29). 



C. Operators on multivectors 

1. The inverse operation 

Given a general three-space multivector 

M = a + w + L-w + Lh (33) 
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we seek an inverse multivector, so that the division operation can be achieved. Firstly, we 
find that 

Z = MM^ = (a + v + LW + Lb)(a - V - LW + ib) = - t;^ + -u;^ - 6^ + 2i(a6 - v • w). (34) 

However this produces a number with a scalar and trivector components, we therefore need 
to multiply by its conjugate in order to create a scalar, giving 

d = ZZ = {a^ -v^ + - b^ f + 4(a6 - v • w)^. (35) 

Therefore we find the inverse 

M-^ = M^Z/d (36) 
= (a — V — tw + t6) (a^ — + — b'^ — 2i{ab — \ ■ w) / d. 

This will fail to produce an inverse if o? = 0, that by inspection requires both 

a'^ + = v'^ + &2 ab^v-w. (37) 

2. The square root 

To find the square root of a multivector M, we need to find a multivector N, such that 

N'^ = M (38) 

where 

M = a + \ + i^ + ib (39) 
is a general multivector. In order to streamline calculations we replace 

M = (a + t6)(l + m + in) (40) 

where by inspection m = "'^^^^ , n = "'^"'''^ and Hq = + b^. Clearly this result is not valid 
if ho = 0, which would imply that a — b — implying a pure vector and bivector, which can 
be found to have the solution sinh^ = — ^±-i/% + l, where go = "^^"'^ and po — • w. 
A subclass of solutions is also produced if po = 0, giving c — d — The case with go = 
in fact appears to have no solution, that is the multivector M — u-\- tU"*-, has no square 
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root, although it can be approximated arbitrarily closely with = e + ^ + ^ + e, giving 
M = u + iU"*- + 2ie^. So for cases with /ig 7^ we seek a multivector such that 

= (c + x + iy + irf)(c + x + iy + ic?) (41) 
= (c^ + - - c/^) + 2(cx - dy) + 2t(cy + cix) + 2t(cci + x ■ y) 
= 1 + m + in. 



We can see that now we have constructed a simplified problem, of the square root of a 
multivector with 6 = 0, and a = 1, but from which we can reconstruct the full square root. 
Inspecting the six linear equations formed from the vector and trivector components we can 
see that we require 

dn + cm cn — dm , , 

where h = + d"^. For the case h = 0, which implies c = d = 0, also implies N = a + tb, 
which has a known square root as for a complex number, calculated later in the appendix, 
and so for more general cases we can assume h > 0. We now need to just find c and d from 



the two remaining equations in Eq. (41) of 



9 ,9 [m^ — n^){c^ — d^) Acdmn 

c -d^ + + —rrr- = 1 43 

cdirn^ — n'^) + (c^ — d'^)m ■ n 



and after multiplying through by h"^ and making the replacement q = "^^^ "'^ and p = n, 
we find 

(c^ - d^) {c^ + d^f + q{c^ -d^) + 2cdp = {c^ + dy (44) 
2cd{c^ + d^ f - 2cdq + (c^ - d^)p = 0. 



We have difficult simultaneous polynomials in c and d, and so it is now convenient to make 
the substitution 

6 6 

c = r cosh - , d = r sinh - , (45) 
2 2 ^ ^ 
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finding 



cosh^ 6 + q + p sinh 6 
sinh 6 cosh^ — sinh 9q + p 



= cosh^ 6 
= 0. 



(46) 



From the second equation in Eq. (46), we find 



1 

^ / g-pcosechg \ * 
V cosh^ 6' / 



We can see that either sign for r satisfies the two equations of Eq. (46) although this will 
only fiip the overall sign of the square root, as expected, giving us a square root for the 
multivector of either sign. The term under the fourth root, q — pcosech.6 in fact remains 
non-negative, and hence r remains real. On substitution back into the first equation in 



Eq. (46) we find the trigonometric equation, 

2q — pcosech6'(l — sinh^ 9) = cosh 9^/q — p cosech^ (48) 
and substituting x = sinh 6 we find the quartic equation in x 

{q - p^)x'^ - p{l + 4g)x^ + [2p^ - q{Aq - 1)) + p(4g - l)x - = (49) 
which has a solution 

X = ^<^+;^+;'^" . (50) 

4(g-r) 

where s = + v/(4g^TFTl6p2 and u = V2^/l6q^ + 12qp^ + p'^ - 4g2 + s(p2 + 4g2)_ ^e 
can see that s remains real, but if we allow a negative square root of s inside the equation 
for u, then we will produce an imaginary term and hence we need to maintain a positive 



square root for s which leaves just the two real solutions to the quartic as shown in Eq. (50). 
Hence we have two distinct 6 = arcsinhx, which imply two distinct square roots, ignoring 
signs. We can see that there is a special case q = p^ that will reduce Eq. (49) to a cubic. 



p{l + 4:p^)x^ + p\4:p^ - 3)x^ - p(4/ - l)x + p = 0, (51) 

which has the single real solution x = —p, giving a single root with 6 = — arcsinhp. If 
p = 0, then the cubic fails and we find two cases dependent on q, with firstly cosh 6* = 2y/q, 



requiring q > \ and secondly 6 = with r = c, where c = y ~^^2 ^^^^^ for g < |. 
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The square root of a + tb is easily found to he f + tg = ±^{\/a + y/h^ + L\/—a + y/ho) 



Hence for M given by Eq. (39), we find the square root 

= ±(c + X + 6y + id) (/ + ig) , 



(52) 



where c, d given by Eq. (45 ) and x, y given by Eq. (42 ). Because multivectors are associative 



we can now find the rational powers M'^/^*, where k,t are integers, which can approximate 
to any fractional power. We find that the square root generally exists, though in some 
cases there is a double root (ignoring signs) and for the special case mentioned we can only 
approximate the root arbitrarily closely. 
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